Thermal contraction in silicon nanowires at low temperatures 
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Abstract 

The thermal expansion effect of silicon nanowires (SiNW) in [100], [110] and [111] directions with 
different sizes is theoretically investigated. At low temperatures, all SiNW studied exhibit thermal 
contraction effect due to the lowest energy of the bending vibration mode which has negative effect 
on the coefficient of thermal expansion (CTE). The CTE in [110] direction is distinctly larger than 
the other two growth directions because of the anisotropy of the bending mode in SiNW. Our study 
reveals that CTE decreases with an increase of the structure ratio 7 = length / 'diameter , and is 
negative in whole temperature range with 7 = 1.3. 
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I. INTRODUCTION 
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Silicon nanowires (SiNW) can be well produced in different directions with different 
l|-|6| It is a promising candidate for nano-device, and has attracted a lot of 
19] The SiNW can serve as basic building blocks for electrically based sensors- 



diameters, 
interest. (7 



field effect transistor 
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16| and logic gates [17|, photovoltaic device [18|], and functional 
networks[19]. The possible stable configurations of the SiNW have been investigated with 
density-functional tight-binding simulations by Zhang et al. in Refs. 
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They found that 



the stability of SiNW is determined by the competition between the minimization of the 
surface energy and the minimization of the surface-to- volume ratio. As one of the important 
thermal properties in SiNW, the thermal expansion plays an important role for these appli- 
cations in the electronic or thermal devices. The thermal expansion effect can be studied in 
the classical molecular dynamics. 21| 

In this paper, we first analyze various phonon vibrational modes in SiNW from the 
Tersoff potential 22] implemented in the "General Utility Lattice Program" (GULP). 23] 
Some important features of phonon modes in SiNW of different directions and with different 
structure ratio 7 = length / diameter are discussed. We then investigate the coefficient 
of thermal expansion (CTE) of SiNW in [100], [110] and [111] growth directions by the 
nonequilibrium Green's function approach [24|], which is a quantum mechanical method and 
automatically includes contribution of all phonon modes. We find that at low temperatures, 
all SiNW studied have thermal contraction effect and the CTE in SiNW [110] is the largest 
one among the three directions. The CTE decreases rapidly with increasing structure ratio 
7, and is negative in whole temperature range with 7 = 1.3. 



II. RESULTS AND DISCUSSION 



Fig. [T] demonstrates that a SiNW with length and diameter (L, D) can be cut from the 
bulk Si crystal by using a virtual cylinder with structure parameters (L, D).[25| We can cut 
SiNW in different growth directions: [100], [110] and [111], by controlling the direction of 



the virtual cylinder. The SiNW disp 
-D) = (3, 1) nm. As discussed in Refs. 



ayed in the figure is in [100] growth direction with (L, 
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it is not important to include H-passivation on the 



surface for thermal property of SiNW. We adopt this approximation which may introduce 
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some small deviation. To depress this small deviation, we do not use these atoms on the 
surface and boundaries. It seems that Si atoms in the figure do not distribute uniformly in 
the figure. That is because the structure is viewed from the upper front direction, which 
is better for illustrating a three-dimensional structure. Other configuration figures in this 
paper may also have this visual effect. 

The frequencies and eigen vectors of phonon modes are obtained from the Tersoff in- 
teratomic potential for silicon implemented in the GULP. The Tersoff potential can well 
describe bonding in covalently bonded solids such as silicon and carbon system, and also 
includes the nonlinear interaction. Its accuracy has been confirmed by comparison with re- 
sults from ab initio density functional theory [261] In the investigation of CTE in the SiNW, 
we will use the clamped boundary condition, where one end (left, for example) is fixed and 



another end (right) is free. For consistency, we a. 
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so apply this boundary condition in the 
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study of phonon, which is different from Refs. 

Fig. [2] shows the vibrational morphology of three important phonon modes for CTE in 
SiNW. In the figure the SiNW is in [111] growth direction with (L, D)—(3, 1) nm. The 
left-most 0.5 nm region (blue online) is fixed (clamped). Arrows (red online) attached onto 
each atom display the vibrational displacement for the atom. [29] The length of the arrow is 
not the exact value of vibrational displacement to show this displacement more clearly. It 
has been enlarged by a constant factor to all atoms. Panels (a) and (b) are the first and 
second order bending modes, which is a lateral mode with all atoms vibrate in the direction 
perpendicular to the axial direction. It is doubly degenerate due to two independent lateral 
vibrational directions. As the SiNW is clamped on the left, the vibrational displacement 
increases from left to right in both modes. In panel (a), the direction of the vibrational 
displacement does not change, while this direction changes in the second order mode in panel 
(b). The bending mode is actually a characteristic property for rod-like systems (eg. SiNW) 
or two-dimensional plane sheet (eg. graphene). Using its relation to the Young's modulus of 
the material, the bending mode has been used to detected the value of the Young's modulus 
of the low- dimensional systems such as single-walled carbon nanotubes 30j or graphene. 311 ] 
For longer SiNW, this mode will gradually turn to the standard flexure mode in a rod- 
like systems. Panels (c) and (d) are the first and second order twisting modes. Because 
of its special rotating morphology, the vibrational displacement increases from inner to 
the surface. This is consistent with the vanishing surface stress condition in a rod. 27J In 
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panel (c), the direction of the vibrational displacement is the same for all atoms, while 
the direction changes for the second order mode shown in panel (d). The twisting mode 
reflects the rod-like structure of the SiNW and is determined by the shear stiffness of the 
SiNW. It can serve to sense the tiny torsion movement of a nano device. This mode can 
be used to extract the internal information of the SiNW, as it carries the information of 
the rigid rotational invariance symmetry. 33M35| Panels (e) and (f) are the first and second 
order longitudinal vibrational phonon modes. In very long SiNW, this mode is actually the 
well-known longitudinal acoustic mode. 

Fig. E] shows frequencies of the three lowest-energy phonon modes of SiNW in different 
directions with different sizes. Frequencies for all phonon modes decrease with increasing 
structure ratio 7 = length / 'diameter . This result indicates that longer SiNW are more 
flexible and elastic. As a result, SiNW with large 7 may be thermally more unstable, and 



possibly show a coiling-up behavior like a single-walled carbon nanotube. [36| Typically, we 
note that the frequency of bending mode in (a) is the lowest one compared with the other 
two modes shown in panels (b) and (c). This is because the bending mode turns to be the 
elastic flexure mode in the limit of 7 — > 00. An important feature of the flexure mode is its 
quadratic phonon spectrum: u cx k 2 . While the other acoustic phonon modes are linearly 
dependent on the wave vector. So the flexure mode (also bending mode) always has the 
lowest frequency. This figure also shows that the bending mode has lowest frequency in 
[100] direction and highest frequency in [110] direction with the same 7. So the SiNW in 
[100] direction are easiest to be bend, which may be due to its regular surface structure. 

We calculate the CTE along the axial direction of the SiNW using the nonequilibrium 
Green's function method, which includes all phonon modes automatically and can be applied 
to very low temperatures as it is a quantum approach. [24] In this method, firstly, the average 
vibrational displacement of atom j, (uj), is calculated through: (uj) = ih~Gj, where Gj is 
the one-point Green's function and can be calculated from its Feynman diagram expansion 
in terms of the nonlinear interaction. We consider the nonlinear interaction as, 

H n = -^ U l U rnU n , (1) 

Ivan 

where the nonlinear force constant is extracted from the Tersoff potential in GULP by finite 
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difference method. Then the CTE can be calculated through its definition, 

d(uj) 1 

where Xj is the lattice position of atom j along the axial direction. We can get the final 
value of CTE for SiNW by averaging oij over all atoms. 

In Fig. HI we compare the CTE in three different directions. All three SiNW have the 
same structure parameters (L, D) = (1, 4) nm. It shows that the CTE in all directions 
are negative in very low temperature region. Then it increases with further increasing 
temperature and reach a saturate value in high temperature limit. The saturate value is 
positive for these SiNW shown in the figure. We mention that before reaching the saturate 
value, all three curves still show obvious increasing behavior at about T = 500 K, which 
implies that the Debye temperature in SiNW is roughly about 500 K. From this figure, we 
also learn that the CTE in [110] direction is distinctly larger than the other two directions. 
While the CTE in [100] direction is close yet a little smaller than that in [111]. 

The above behaviors for the CTE at different temperatures and in different directions 
can be understood from the various phonon modes in SiNW discussed in the previous part. 
Because of different vibrational morphology, these phonon modes have very different con- 
tribution to the CTE. The bending mode has negative effect on the CTE of the rod-like 



SiNW. 
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38] The twisting mode only takes up some degrees of freedom, without making 
any contribution to CTE. The longitudinal vibration mode has positive effect on the CTE. 
So the CTE is mainly determined by the competition between the bending mode and lon- 
gitudinal mode. Since the bending mode has the lowest energy among all phonon modes, it 
is the only mode excited at low temperature, leading to negative CTE. With the increase 
of temperature T, the longitudinal mode will be excited gradually and competes with bend- 
ing mode. This competition slows down the decreasing of CTE, and results in a minimum 
point. As the longitudinal mode becomes more and more important, CTE will increase 
with further increasing temperature. In high temperature limit, CTE reaches a saturate 
value after sufficient competition between all phonon modes. The positive saturate value 
in these SiNW implies the leading contribution of longitudinal mode for these SiNW. Also 
the bending mode is anisotropic with the frequencies in [110], [111] and [100] directions as 
53.1, 41.2, 39.3 cm -1 , respectively. This anisotropic property leads to the anisotropic effect 
of CTE, i.e. a[i 10 ] > a [ul] > apo]- 



We now turn to the discussion of the size dependence for CTE in SiNW [111]. There 
are similar effects for the other two directions. Fig. [5] shows the dependence of CTE on 
the structure ratio 7 in whole temperature range. The value of CTE decreases quickly with 
increasing 7. This is because the frequency of bending mode decreases quickly with increas- 
ing 7 (see Fig. [3]). As a result, the negative effect from bending mode on CTE increases. 
For SiNW with 7 = 0.25, the transition of CTE from negative to positive happens around 
120 K, and CTE reaches a saturate value of 10 x 10~ 6 K _1 in high temperature limit. A 
large positive saturate value of CTE indicates that the longitudinal mode makes dominant 
contribution rather than the bending mode in this SiNW. The transition temperature in- 
creases from 120 K to 270 K with the change of 7 from 0.25 to 1.0. The saturate value of 
CTE for SiNW with 7 = 1.0 decreases for more than 50% compared with that of 7 = 0.25. 
This result is reasonable and has also been observed in single-walled carbon nanotubes, 24j 



which has similar one- dimensional rod-like structure. It seems to be a common phenomenon 
in one- dimensional system that the CTE is much smaller in longer and thinner structures. 
When 7 = 1.3, CTE is negative in whole temperature range, which indicates that the SiNW 
with these structure ratio undergoes thermal contraction at all temperatures. For these type 
of SiNW, the bending mode has played an dominant role over the longitudinal mode. 



III. CONCLUSION 



To summarize, we have studied the thermal expansion of SiNW in different growth di- 
rections by the nonequilibrium Green's function method. We find that at low temperatures, 
all SiNW studied exhibit thermal contraction effect due to the lowest energy of the bending 
mode which has negative effect on CTE. The CTE in SiNW [110] is larger than the other two 
directions, because of higher frequency of bending mode in SiNW [110]. The CTE decreases 
quickly with increasing of the structure ratio 7, and is negative at all temperatures with 
7 = 1.3. 
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FIG. 1: (Color online) Use a virtual cylinder to cut a SiNW from bulk Silicon. 
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(a) 



(c) 




(e) 




(d) 
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FIG. 2: (Color online) The morphology of different phonon modes in SiNW. (a)/(b): the 
first/second order bending mode; (c)/(d): the first/second order twisting mode; (e)/(f): the 
first/second order longitudinal mode. 
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FIG. 3: (Color online) The frequency of three lowest-energy phonon modes in SiNW with different 
size. (a), bending mode, (b). twisting mode, (c). longitudial mode. 
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FIG. 4: (Color online) CTE of SiNW in three directions [100], [110], and [111]. 
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FIG. 5: (Color online) CTE depends on T, in [111] SiNW with different structure ratio 7 = 
length / diameter . 
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